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Abstract 

We consider the resonant production of fermions from an oscillating axial 
background. The classical evolution of the axial field is given by that of 
a massive pseudovector field, as suggested by the renormalizability of the 
theory. We look upon both the massive and the massless fermion production 
from a perturbative point of view. We obtain the corresponding spectrum and 
angular distributions for the different spins or helicities in the particular case 
of a spatial-like axial field. We also extend our study to the non-perturbative 
regime in the massless case and compare the results with the perturbative 
ones. 



1 Introduction 

The parametric resonance phenomenon plays a fundamental role in the modern 
theories of reheating after inflation [I]. In such theories, after a period of inflation, 
driven generally by the potential energy of the slowly rolling inflaton field, the 
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inflaton experiences a process of coherent oscillations around the minimum of the 
potential. Due to its coupling to the rest of the fields, those oscillations can give 
rise to the production of exponentially large amount of bosons and subsequently 
to fermions through the decay of newly formed bosons, whose energy spectrum is 
characterized by resonance bands. Although, it is an efficient mechanism to convert 
the inflaton energy to the intermediate bosons, it is still far from complete and we 
still require to invoke the old idea of reheating where the bosons decay to various 
fermions and that subsequently leads to the thermalization of the universe. The 
standard picture of preheating has not been tested for fields other than inflatons, 
and therefore it is still unexplored what would be the effect of other fields that could 
be present in the early universe such as dilaton, axial or moduli predicted by string 
theory. It is a general belief that there would still be an amplification to one of 
the dominant fields depending upon the values of coupling constants but it has also 
been shown that for two scalar fields the system becomes chaotic ||. 

The production of particles from classical sources has other important conse- 
quences apart from reheating. In fact, the same phenomenon is responsible for the 
amplification of vacuum fluctuations in the generation of gravitational waves || . It 
has also been proposed as a mechanism for the generation of primordial magnetic 
fields in the context of string cosmology Q. In this latter model, a dynamical dilaton 
field plays the role of the external source. In addition, many other examples of par- 
ticle production based on alternative gravity theories can be found in the literature, 
such as Brans-Dicke ||, higher dimensional models ||, etc. 

In this paper we do not intend to give a holistic picture of preheating; rather 
we explore a different mechanism for the amplification of vacuum fluctuations, it 
is based in the presence of an oscillating background axial field. Axial (or anti- 
symmetric) fields appear naturally in the bosonic sector of the low-energy string 
effective action together with graviton and dilaton fields [0, |§. m f&ct, recently 
some explicit solutions with non-trivial axion fields have been obtained in the con- 
text of homogeneous, isotropic and anisotropic models in string cosmology ||, [10| . 
These antisymmetric fields can be interpreted as torsion in pseudo-Riemannian ge- 
ometry || and, as we will show, this fact provides a natural way of coupling them 
to other matter field. This is also the case of supergravity theories which contain a 
non-symmetric part in the spin-connection that is determined by the gravitino field 
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The presence of general chiral fields can give rise to a different kind of production 
from the usual scalar fields. In particular, the breaking of parity invariance could 
lead to a different production of left and right fermions. In the present case, this will 



2 



not happen because our model is parity invariant. In addition, as far as rotational 
invariance is broken by the axial field, we will obtain an anisotropic production of 
particles. In the usual picture of reheating it is expected that such initial anisotropy 
is rapidly erased by subsequent decays of the fermions produced in arbitrary direc- 
tions. On the other hand, as far as the background field could have a non-vanishing 
angular momentum, we should also pay attention to the spin distribution of the 
produced particles. Again in our case, we will consider a very simple homogeneous 
background with vanishing total angular momentum and we will not generate any 
net spin. 

The paper is organized as follows: in Section 2, we study the classical equations 
of motion for the axial field and its coupling to fermions. Section 3 is devoted to the 
perturbative calculation valid in the limit of small axial field for both massive and 
massless fermions. In Section 4 we extend our calculation to the non-perturbative 
regime in the masless limit and compare our results with the previous ones. In 
Section 5 we include the main conclusions of the work. 



2 Minimal coupling and the axial field dynamics 

The torsion field can be seen to be minimally coupled to fermions by means of 
the Einstein Equivalence Principle [p~2| , |T3f . Thus, the Dirac lagrangian for a single 
fermion field wih mass m in the presence of curvature and torsion is given by: 

£ = (iy + fi M + igStfs) -m)iP (1) 

where S p = e pv \ p T^ uX is the torsion pseudotrace, Q p the spin-connection and the 
coupling constant is fixed g = 1/8, (we will work with a general g in order to extend 
our results to other models in which S p is not interpreted as torsion). Despite the 
presence of the 75, the lagrangian is parity invariant since S p is a pseudovector. As a 
consequence of this coupling, the axial field can act as a source of fermion creation. 
Notice that this coupling is dictated by the Einstein Equivalence Principle, unlike 
the usual inflaton couplings which are introduced ad hoc. The problem of fermion 
production has received attention only very recently in the works of Greene and 
Kofman |TJJ and Baacke et al. |T5|] . They have studied for the first time the 



problem of reheating of fermions from the coherent oscillations of the inflaton field. 
They showed how the limit on occupation number imposed by Pauli blocking was 
saturated by the resonant production. 
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The first problem we have to deal with regarding the creation of particles is that 
of the classical dynamics of the torsion field. There are several models proposed 
in the literature for which torsion appears in the classical action only as a mass 
term, thus for instance: the Einstein-Cartan theory considers the Einstein-Hilbert 
action but replaces the scalar curvature built out of the Levi-Civita connection by 
the scalar curvature built out of an arbitrary connection with torsion. In this theory 
torsion is a non-propagating field that appears in the lagrangian simply in a mass 
term of the form MpS 2 , with Mp the Planck mass. This is also the case of the 
bosonic sector of the low-energy string effective action. However, if we consider the 
coupling of torsion to quantum fermionic fields in ([!]), the conditions of unitarity 
and renormalizability automatically requires the existence of a kinetic term for the 
axial field, together with the previous mass term, i.e, they impose the axial field to 



behave like an abelian massive axial-vector [16, 17]. In fact, the vacuum divergences 



generated by the above lagrangian in are given by: 

sus] = (Aj/ A (-^ + ^) (2) 

where = d^Sy — d^S^, A = N e + \og(fi 2 /m 2 ) with the poles parametrized as usual 
in dimensional regularization by N e = 2/e + log47r — 7 and fi is the renormalization 
scale. Notice that there is no S 4 term although it has the same dimension. When 
instead of a single fermion we have more than one, the previous action is replaced 
by the sum of the actions corresponding to each single fermion. Therefore by an 
appropriate renormalization of mass and wave function we can take as the classical 
lagrangian for torsion that of the abelian massive gauge field. This is the minimal 
lagrangian for torsion that ensures the renormalizability of the fermionic sector. This 
model is analogous to the inflaton model with a quadratic potential, but replacing 
the scalar field by a pseudo-vector. Notice that unlike the vector abelian case, the 
presence of massive fermions is incompatible with a gauge invariance associated to 
the axial field. The corresponding equation of motion is nothing but the Proca 
equation, i.e: 

d^ v + m 2 s S u = (3) 

with m s the mass of torsion. It is now easy to find solutions for this equation. 
In particular, when torsion is spatial-like and only depends on time we have the 
following oscillatory solutions: S l (i) = S l sm(m s (t — to)) with S l constants. As is 
well-known, this kind of periodic functions can give rise to the parametric resonance 
pheonomenon. 
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As in the inflaton models, the production of particles induces the damping in 
the torsion oscillations. This can be described either by including a friction term in 
the equations of motion or by means of the vacuum polarization correction to the 
torsion mass. The vacuum polarization correction has been calculated in a previous 



paper 18] and from its imaginary part we get that: 



ImII(m s ) q 2 m 2 , , 

oc (4) 



m q m. 



where m? is the mass of the fermion squared or in the case of several fermions 
present, it is given by m 2 = Yji m 1- Accordingly the classical evolution can be taken 
as: 

S*(t) = & sin(m s (t - t ))e~5 r (*-*°) (5) 

For the sake of simplicity we do not consider the expansion of the universe and 
accordingly all the calculations will be done in Minkowski space-time. We will 
also simplify further and take S l pointing in the z-direction. As usual in particle 
production calculations, we will impose our initial conditions on the background 
fields in such a way that it is possible to define asymptotic in and out vacuum 
states. With that purpose we will assume in next section that when t — ► — oo the 
axial field is zero, i.e we will take S l = for t < t and S 3 (t) given by (0) for t > t , 
as we will show the result will not depend on the particular value of t but only on 
the amplitiude of the oscillations S l and decay rate T. 



3 Perturbative approach 

Massive fermions 

Let us consider the Dirac equation in the presence of the axial field: 

(17^-5^75 -m)* = (6) 

It is to be noted that the form of the coupling does not allow us to disentangle the 
different spin modes, unlike in the case of Greene and Kofman. In addition, it cannot 
be reduced in an easy way to the known form of Mathieu equation or Lame equation, 
which has the known stability and the instability bands. During instability bands 
the occupation number grows exponentially fast in the case of bosons. This keep us 
from finding the exact solution of the equation. However, it is possible to perform a 



5 



perturbative analysis which provides part of the resonant behaviour of the equation. 
There exists various approximation schemes to apply for the bosons as well as the 
fermions; we will follow that in [19|. Multiplying (|6|) with (ij^d^ + gSf/y^^ + m), 



we get 

+ i^-fjsd^Su + igSui-f, - g 2 S 2 + 2gm S l5 + m 2 } * = (7) 

Since there is no dependence on position in the Dirac operator, we can write the 
general solution \I/ as: 

*(t 1 2)=^,.(t)e 4fe (8) 

where s is a spin index. As we mentioned before we will take = (0,0,0, S3). 
Eq.(0) reduces to: 

(d 2 + k 2 + ^47°7 3 75 - 2gk 3 S 3l5 + m 2 + g 2 S 2 3 + 2gmS 3l 3 l5 ) ^ s (t) = (9) 

We introduce the plane- wave solutions for massive spinors with positive and negative 
frequency, and with a definite spin projection along the z-direction: 

U^ s (t,x) = -t=u(k,s)e^-^ (10) 



Vr(t,x) = -=v(k,s)e- iks+iUit (11) 
2oj 



where s — ± and uj 2 = k 2 + m 2 , u(k, s) and v(k, s) satisfy the following normalization 
conditions: 

u ] (k,r)u(k,s) = 2u8 rs (12) 

v' (k,r)v(k, s) = 2u5 rs (13) 

u j (-k,r)v(k,s) = (14) 

v\-k,r)u(k,s) = (15) 

where u(k, s) and v(k, s) are defined as: 

u(k, s) = \/uj + m ^ ff ^ s ^ j , v(k, s) = \fuT+m y ^+m^ s ^ (16) 
and the Dirac spinors have the simple form: 



We will work in this section in the Dirac representation for which the 7 matrices 
appear as: 

U- :?)■*- (A o 1 ) <»> 

Let us consider a given solution of the Dirac equation ip m that initially, i.e. t — > 
—00, behaves like a plane wave with positive energy and spin s. In the remote 
future, where the interaction is switched off, such a solution will evolve into a linear 
combination of positive and negative frequences with a certain probability for spin 
flip, that is: 



f ) ^ E «Sm,0&/ + Pi~V-i« ( 19 ) 



/ — X 

s/ 

The Bugolubov coefficients for fermions satisfy the following relation: 

E(KJ 2 + l/%J 2 ) = i (20) 

si 

It is possible to transform the differential equation into an integral equation that 
allows a perturbative expansion: 

€(M>^ + -f M^in(u(t-t')W£(t',x)dt' (21) 

Cl> J —00 

where the interaction term is given by: 

Mg = g^Sf - 2p fc 3 g 3 7 5 + fom gys + + igS^Vjs (22) 

7 II 7/7 71/ V 

From (P^D , we get for the Bugolubov coefficients to first order in perturbation theory 
the following expression: 

, j roc 

, = ,ibf = / v ] r Mrur e~ i2uJt dt' (23) 

t^kss' _ks' Y ks 4^2 J_ oo -ks' k ks \ I 

where we have substituted ■0i n in the r.h.s of (|2l|) by the lowest order solution U% a 
and we have made use of the limit t — > 00 of equation (pT]). The corresponding 
number of particles in a given spin state is given by the sum: 

^ = EI/%J 2 (24) 



7 



Now we discuss the contribution of each term present in the potential. First we notice 
that the diagonal terms / and IV do not contribute to first order in perturbation 
theory. The V term gives rise to: 

^^tVAe. = 2gSs(h T iki) (25) 

where — signs applies for s' — — and s = + and the + sign for opposite spin 
directions. For s' = s the term vanishes. The II term gives rise to 

u + m- ^ + m j (26) 

Finally the /// term yields: 

2gmv i _ ] . s S 3 j 3 j 5 uj; s = -igmS 3 k 3 (27) 

Notice that the contributions from // and III cancel each other and finally only the 
term V contributes to the Bogolyubov coefficients. Performing the Fourier transform 
implicit in (^3f), and using the spherical coordinates in momentum space (k, 6, </>), we 
can write the total number of particles produced with spin s = ± and momentum 
k as: 

p 

N u = N i _=g 2 \S 3 {Lo)\ 2 -^6 (28) 

Uu 

with S 3 (lu) the Fourier transform of S 3 (t) and 6 the angle between k and 5*. As 
commented before, the number of produced fermions with spin up equals the number 
of fermions with spin down and as a consequence there is no net spin creation. The 
total number of fermions is just twice the above result. 

We see that the production is not isotropic, with a maximum in the orthogonal 
directions to S (6 = vr/2) and no forward-backward production (9 = 0,tt). In 
Fig.l we have plotted the spectrum in the lower curve and we observe the resonance 
close to k = m s /2 which corresponds to the decay of the axial field quanta into 
fermion-antifermion pairs. 
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Massless fermions 

As we have seen, it is not possible to separate out the positive and negative spin 
states in the massive case, but in the massless case it is possible to do that by choos- 
ing a different representation for the Dirac matrices, the Weyl representation. Since 
the solutions of the Weyl equation in the absence of external sources can be written 
as eigenstates of the helicity operator, it is useful to consider those eigenstates for 
the positive and negative frequencies solutions. We shall derive the total number of 
massless fermion production and we shall see that the summation of both chiralities 
gives rise to total number of fermions and it is equal to the one obtained in the 
previous section taking m — > 0. The form of the plane-wave solutions with positive 



and negative frequency is that given in (fTT[) , but now the form of the spinors is given 
by: 

= ^ ( ^ V H- = ^ ( „ %, ) (29) 



^^(T)'^^(^)) (30) 

where 

a + (k) = -6 + ffl = ( ™j« ) , a4k) = -MS) = ( - Z}f* ) <»> 

are eigenstates of the helicity operator with eigenvalues +1 and —1 respectively. 
The normalization conditions are the same as the massive ones. In the Weyl repre- 
sentation the gamma matrices are: 

v-U oO'^U'o')— (i-i) < 32 > 

In the present case, the initial positive frequency solution will evolve into a linear 
combination of positive and negative frequency modes, but since the perturbation is 
diagonal in chirality space there is no mixing between positive and negative chirality 
modes, i.e: 



+ (33) 
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Figure 1. -Number of particles against |fc| for gS = 1/8 m s and T = 0.2 m s . For 
massive fermions m = 0.3m s 

In the massless case, the contributions in the III and IV terms are not present and 
again the diagonal terms do not contribute. As a consequence only the V term is 
relevant. Hence the total number of fermions with momentum k and helicity s is 
given by the following expression: 

Ni = \P i f = fsin 2 0\S 3 (u>)\ 2 (34) 

The total number with both the spins would be just twice the above expression, 
which is the same as the massive case. In Fig. 1 we have represented the spectrum 
for massless fields in the upper curve. It is to be noted that the production of 
massive fermions is, as expected, suppressed by the mass term. We also observe the 
peak at k — m s /2. 
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4 Non-perturbative results 



In the previous section we have obtained the spectrum and angular distribution of 
the produced particles up to first order in perturbation theory. Higher order terms 
are expected to be suppressed and only will give rise to smaller peaks at higher 
energies. However, the non-perturbative effects can have a more important effect on 
the results. As shown in the scalar case, the features of production essentially deviate 
from the pertubative result. Since in our case the equations cannot be reduced to 
any known form, we will solve them numerically. For the sake of simplicity we 
will only consider the massless case in which it is possible to disentangle the two 
spin states. We will also assume that the behaviour of the background field is the 
following: for t < and t > nT S3 = S a constant, with n an integer and T = 27t/m s 
the period of the oscillations. For < t < nT S3 = S cos(m s t). It has been checked 
that the final spectra do not depend on n. In order to reduce @ to a second order 



form we make the following ansatz on the solutions [14 



i> = {ifi + g$l,)h s {t)J kS W s ± (35) 
where for positive (negative) + (-) energy modes with spin s we take: 



W s+ = *f , W- = ( U (36) 





In this section we will work again in the Weyl representation. With this ansatz, the 
equations of motion can be written as a single equation for the function . 

(d 2 + k 2 - isgS 3 + 2gk 3 S 3 + g'S^f^t) = (37) 

In addition the ip spinors will give rise asympotically (t — > ±00) to eigenstates of the 
the third component of the spin operator. In this sense, they have the appropriate 
form to describe the creation of fermions with definite spin along that direction. The 
normalized spectrum of the total number of particles created with spin s is given by 
(see 0, |2TJ for details): 



k 2 + k 2 1 
% = fr + Jfr +\gSk 3 sm 2 (d 2 ) {lmF ^ {T))2 
Psin 2 (£) (Imi^(T)) 2 



k 2 + g 2 S 2 - 2gSkcos6l - (ReFg s (T)) 5 
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Figure 2. -Normalized number of particles against |fc|. a),b) and c) correspond to 
gS = 0.1,1, 10m s respectively with 9 — tt/2 and d), e) and f) to the same amplitudes 
but for = 7r/4. 
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e) 

Figure 3. -Normalized number of particles against 9. a),b) and c) correspond to gS = 
0.1, 1, 10m s respectively with \k\ = .lm s and d), e) to gS = 10m s with \k\ — 1, 10m s 
respectively. 
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where F^ s is a solution of equation (|37| ) with initial conditions F^ s (0) = 1 and 
Ft (0) = and we have used the fact that cosc^ = HeF^ (T). The numerical com- 
putations show that the number of fermions produced with spin up is the same as 
those with spin down, as expected there is no net spin creation in this particularly 
simple model. We see that sin 2 9 angular dependence that we obtained in the per- 
turbative calculation is already present in the above result for S ^> k. Although it 
is expected to be modified by the factor involving (T). In Figures 2 and 3 we 
have represented the spectra against |fc| and 9 for different values of the amplitudes. 
First we notice that for small amplitudes of the oscillations S, there appears a peak 
at k = 0.5m s in agreement with the perturbative calculation. However we find new 
peaks also for small amplitudes. For large values of the amplitude the departure 
from the perturbative results is apparent. We can also observe the appearance of 
the resonance bands in the spectra whose structure strongly depends on the angle 
9. In addition, the angular dependence is also highly sensitive to the energy k of 
the fermions and the amplitude S of the oscillations. In particular we find that for 
large amplitude of the oscillations, strong peaks appear at some particular values of 
the angles, indicating that the production mainly takes place in those directions. 



5 Conclusions 

In this work we have studied the resonant production of fermions from a classical 
oscillating axial background. The dynamics of the classical background is deter- 
mined by the renormalizabililty of the quantum theory of fermions coupled to axial 
fields and is given by the Proca equation. The oscillatory solutions to this equation 
are obtained and the corresponding phenomenon of resonant particle production is 
analyzed both from a perturbative and a non-perturbative point of view. In the 
perturbative case, we have considered both the massive and massless fermion cases 
for a simple background in which the axial field is oscillating in the z-direction, and 
we have obtained the spectra and angular distributions of the produced particles. In 
the non-perturbative case, we consider only the massless case and study the equa- 
tion numerically. We have checked that the non-perturbative results deviate from 
the perturbative ones with the appearance of the typical resonance bands, that in 
our case depend both on momenta and angle. 

Concerning the possibility that this mechanism could have been relevant in the 
early universe, we should mention that the generation of an initial anisotropic axial 
field that later on could have evolved according to (0) has been explored in f22 |. 
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There it is shown that for axial fields interacting with photons, there is a possibility 
of spontaneous symmetry breaking along spatial directions for the axial field. As a 
consequence, that field would acquire some spatial-like vacuum expectation value. 
In addition, as we mentioned in the introduction, axial fields are present in the 
context of string cosmology and therefore it would be very interesting to study if 
the resonant phenomenon studied in this work could affect the standard picture of 
reheating after inflation. 
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